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We investigate the influence of diffraction on the statistics of energy levels in quantum 

systems with a chaotic classical limit. By applying the geometrical theory of diffraction we 
show that diffraction on singularities of the potential can lead to modifications in semiclassical 
approximations for spectral statistics that persist in the semiclassical limit h — > 0. This result is 
obtained by deriving a classical sum rule for trajectories that connect two points in coordinate 
space. 
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1 Introduction 

Correlations in the spectra of quantum systems with a chaotic classical limit are expected to 
follow random matrix theory. This is the content of the random matrix hypothesis [[l], [2| , which 
has found confirmation by numerical investigations in many systems. Theoretical support for 
this conjecture has been obtained by semiclassical approximations based on periodic orbits |3|. [| 
as well as field theoretical methods ||. In particular, it has been shown by the semiclassical 
method that the leading asymptotic behaviour of two-point correlation functions for long-range 
correlations agrees with results of random matrix theory. This result has been obtained by 
using mean properties of periodic orbits as expressed by the sum rule of Hannay and Ozorio de 
Almeida || . It is expected that a proof of the random matrix hypothesis would be possible if 
certain correlations between periodic orbits were known |7[j. 

The semiclassical analysis is based on the trace formula for the density of states ||. This is 
a leading order approximation, as H — > 0, in terms of the periodic orbits of the corresponding 
classical system. For most systems this approximation is not exact, and there are corrections in 
higher order of h. A particular important correction occurs if the potential has discontinuities 
or singularities which lead to wave diffraction. This is the case in many standard examples of 
chaotic systems. It leads to additional contributions to the trace formula in terms of so-called 
creeping orbits or in terms of trajectories that are closed or connect two points in coordinate 
space |9|, [H], [U], . In the present article we investigate systems in which the diffraction occurs 
at point-like objects, for example at a delta-like singularity of the potential, at a magnetic flux 
line in a two-dimensional system, or at a corner in a billiard system. For these systems we 
examine whether diffractive corrections to the trace formula can have an influence on spectral 
statistics that persists in the semiclassical limit h — > 0. The analysis is done by deriving a 
classical sum rule for the orbits that arise in the geometrical theory of diffraction. With this 
input the diffractive corrections to the diagonal approximation for the spectral form factor are 
calculated, and it is shown that these corrections in general do not vanish as h — > 0. 

2 Sum rules for transient orbits 

Semiclassical approximations for the Green function G(q b , q a , E) of a quantum system involve 
classical trajectories that go from q a to q b at energy E. In general there is an infinite number 
of these trajectories. Moreover, in a chaotic system the number of trajectories that connect 
the two points in a time less than T increases exponentially with T. In this section we use 
the ergodic property of chaotic systems in order to obtain a sum rule for orbits connecting two 
points in coordinate space. We call these trajectories transient orbits. 

Consider a particle with energy E that starts at time t = at a point q a in coordinate 
space. The classical probability density for the particle to be found at time t — T at a point 
q b is given by 

P(n . Tn m Id f PaS(E-H(q a ,p a ))6(q(T)-q b ) 
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where the function q(t) denotes the position of the particle with initial position q a and momen- 
tum p a at t = 0, and / is the number of degrees of freedom of the system. The integration in the 
numerator extends over all initial momenta corresponding to energy E, and the denominator 



gives the normalisation constant. The probability density is only different from zero if there 
are classical trajectories that connect the points q a and q b in time T. 

Let us assume that the particle moves in the field of a vector and a scalar potential and 
that the Hamiltonian is given by 



ff(«.P) = ^(p-^A(,)) 2 + y(,) 



(2) 



In the following we express the transition probability density for this Hamiltonian in terms of 
classical trajectories. 

The numerator in ([!]) is evaluated by introducing a local coordinate system in the vicinity 
of a trajectory where one coordinate is along the trajectory and the others are perpendicular 
to it, 
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The sum is over all trajectories from q a to q b at energy E, v a and Vj, are the velocities at q a and 
q b , respectively, T 7 is the time along the trajectory 7, and the prime at the determinant denotes 
that the matrix of derivatives involves only the coordinates orthogonal to the trajectory. 

The denominator is evaluated by changing the integration over momentum into one over ve- 
locity and evaluating it in hyperspherical coordinates. The remaining integral over coordinates 
gives one, and one obtains 
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where S^' = 2tt^ 2 /T(f /2) is the surface area of an /-dimensional hypersphere with unit 
radius. Altogether the result is 
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In an ergodic system the property that typical trajectories fill out phase space uniformly 
can be used to describe the probability density for long times T. In order to apply ergodicity 
we smooth the singular function P over some range of the final coordinate 



(P(q b ; T, q a , E)) E = J ^ q' b P(q' b ; T, q a , E) 5 £ (q b - q' b ) 
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where 5 £ is a smoothed delta-function whose width is parameterised by e, and which is nor- 



malised to one 



/ dfq5 £ (q-q b ) = l 
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It is assumed that the width e is classically small so that the potentials do not change appre- 
ciably within this region. The integral in (|^) results in the replacement of the second delta 
function in the numerator of (|l|) by 5 £ . 



In order to describe the properties of the probability density P we employ the property of 
an ergodic system that the time-average of a function along a trajectory is equal to its phase 
space average for almost all trajectories 

hm - / dt f(q(t),p(t)) = J - — rr , 8 

t^ooT J To Jdtqd}pd(E - H(q,p)) 

and we choose for / the smoothed delta- function S E . The phase space average over this function 
is evaluated similarly as before and yields 
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where S(-E') is the volume of the energy surface in phase space. After integrating @ with @ 
over the initial momenta and dividing by the normalisation constant one obtains the following 
result for the probability density 

hm - / dt (P(q b ; t, q a , E)) e = { b) + 0{e) . (10) 

This property of P has the following interpretation. When the first term on the right-hand 
side is multiplied by a volume element d?q it is the ratio between the volume of the part of the 
energy surface in phase space corresponding the element d^q around q b and the volume of the 
total energy surface. Thus the probability of a particle to be found in a neighbourhood of q b is 
equal to the relative volume of the energy surface of this neighbourhood. 

Equation ( |T0"D gives an expression for the probability density P that is smoothed over an 
^neighbourhood of the final point in an ergodic system. Using (||), it is a sum rule for classical 
trajectories that start at q a and visit the ^-neighbourhood of q b . In order to obtain a sum rule 
for trajectories that hit q b one has to take the limit e — > 0. Here one faces the problem that 
one would like to interchange the two limits e — > and T — > oo, a problem that occurs also in 
the derivation of sum rules for periodic orbits || . 

One can argue that the two limits can be interchanged if the classical trajectories do not have 
conjugate points. Then the terms [det\dq b /dp a )Y 1 decrease exponentially with the transition 
time Tj. This means that the contribution of a single trajectory to the probability density is 
exponentially small for long times. The integral over the probability density j T dt P(q b ; t, q a , E) 
is a discontinuous function of the final coordinate q b , but its variation inside the e-environment 
of q b becomes very small for long times T, and the two limits e — > and T — > oo can be 
interchanged. 

This argumentation does not apply, however, if there are caustics of the classical motion 
inside the ^-neighbourhood. On a caustic the term [det' {dq b / dp a )Y x is divergent. This diver- 
gence is integrable so that the smoothed version of the sum rule (|10"D is still valid. It might be 
possible that in the generic situation where the point q b itself does not lie on a caustic, the two 
limits can still be interchanged. We point out, however, that for applications in semiclassical 
approximations quantum mechanics provides a natural smoothing of the probability density, 
so that only a smoothed version of the sum rule is required. 

In cases where the two limits can be interchanged we obtain the following sum rule for 
trajectories from q a to q b at energy E 
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and the differentiated version is 
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(12) 



where it is implied that the left-hand side has to be smoothed over some small time-interval 
AT in order to obtain a smooth function. 

In the following we give several variants of the sum rule that can be obtained in an analogous 

way. 

• For chaotic area-preserving maps on an (2/)-dimensional unit torus, (q n +x,p n +i) = 
h(q n ,p n ), the corresponding result is 
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where the sum extends over all points (qo,Po) for which q = q a and q n = q b . In appendix 
[A] it is shown that this sum rule is exact for cat maps for all n > 0. 



For billiard systems the sum rule (|12|) can be expressed in a geometrical form. In two 
dimensions the result is 
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where A is the area of the billiard, L 7 is the geometrical length of a trajectory and 
A^i2 = Pa {dq^/dpa ) is an element of the stability matrix, scaled by p a in order to be 
energy independent. Eq. flI3| ) expresses that a particle is equally likely to be found 
anywhere in the billiard if it travels a sufficiently long distance L. Consequently, the 
probability density that it is found at some point q b is equal to A" 1 . 

For applications in the next section one requires sum rules for a subset of transient orbits 
for which the angular orientations of the initial and the final velocities, <f) a and 4> b , are 
fixed, 



E 

0a.0t, fixed 



det' 



dq b 
dp a 



S(T - T 7 



[mv r . 



\f-i 



(mv b ) 



/-i 



X(E) 



dQ n dQh 



T^oo. (15) 



Here the sum extends over all transient trajectories for which the initial and final angular 
directions lie in solid angle elements dfi a and dfl b , respectively, around the directions <p a 
and <p b . 

In integrable systems a similar sum rule can be obtained, if the ergodic average is per- 
formed over an invariant torus instead of the energy surface. 
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where J are the action variables of the system. This equation has the following inter- 
pretation. The left-hand side is a sum over all trajectories for which the initial angular 
direction lies in a solid angle element dQ a around the direction <f) a . This direction (f) a , 
the energy E and the initial position q a determine the action variable I a = I(q a ,p a ) and 
thus the torus on which the motion occurs. The right-hand side is a summation over all 
final momenta p b for which the points (q b , p b ) lie on this torus. 

In order to obtain a sum rule for all transient trajectories one has to integrate (|16D over 
all initial directions. In contrast to the chaotic case this is not always possible, since the 
integral can be divergent due to caustics. In cases where it is possible one has 
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This sum rule can be verified, for example, for a /-dimensional rectangular billiard with 
side lengths <2j, i = 1, . . . , /. There the number of transient trajectories with length in 
an interval dL around L is given asymptotically for large L by <S^ L^ 1 dL/V where 
V = nr=i a i * s the volume of the billiard, and L = v a T. The absolute value of the 
determinant on the left-hand-side of ([[7]) is (L/ijnVa))^ 1 , the actions are U = Oj \pi\/Ti, 
and the sum on the right-hand-side is (2tt)^ /V. Altogether one obtains on both sides 
V -1 . Incidentally the sum rule for chaotic billiard systems yields the same result in this 
case. The reason is that for rectangular billiards the determinant on the right-hand side 
of (HTD is the same for all tori. 



The sum rule ( JE2| ) contains implicit information about the number of transient orbits. Let 
us assume that due to the exponential sensitivity of trajectories on initial conditions one has 
(\det'(dq b /dp a ) | _1 ) ~ d exp(— hT) where d and h are constants, and the averaging is performed 
over trajectories with T 7 p» T. With this assumption one obtains the following asymptotic law 
for the number of transient orbits, 

p(T) = ^1 ~ c exp(hT) , T - oo , (18) 

where N(T) is the number of transient trajectories with T 7 < T, p(T) the density of these 
orbits, and c is a constant that is determined by d . For Riemannian manifolds with no conjugate 
points one can show that limx^ooT -1 logjV(T) = h is the topological entropy of the system 
[|13| . Eq. (|18D with topological entropy h would imply that the number of transient orbits grows 
by an order of T stronger than the number of periodic orbits for which p po {T) ~ exp(hT)/T. 
One arrives at a similar conclusion for the number of transient orbits if one considers systems 
in which a code for these orbits exists JTT], [T3|, |T5|1 . For cat maps the law (|18|) is proved in 



appendix [S|. 

3 Influence of diffractive orbits on spectral statistics 

Diffraction of quantum wave functions on singularities of the potential leads to corrections in 
semiclassical expansions in terms of classical trajectories. In approximations for the density 
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of states there are additional terms besides periodic orbits. These are expressed in terms of 
diffractive orbits. We concentrate on cases where the source of diffraction is point-like for which 
some examples are given below. In these cases the diffractive orbits are trajectories that start 
from and return to the source of diffraction n times, where n is an arbitrary positive integer. 
They are composed of a sequence of arbitrary n transient orbits of the last section for which 
the initial and the final points of the trajectory are identical and are located at the source of 
diffraction. Alternatively, they can be considered as closed trajectories that are scattered n 
times on the source of diffraction. 

Within the geometrical theory of diffraction (GTD) the contribution of all diffractive orbits 
with n scattering events to the density of states is given by || [10], |TT | 
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where 
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The expression (|20| ) is the semiclassical contribution of a transient orbit 7, to the Green function, 
and S 7i is its action and u 7i the number of conjugate points along it. Furthermore, v a and v b 
are the initial and final velocities. They are identical since the initial and final point of the 
trajectory are identical, and we write them in the following without subscript. 

A diffractive orbit is composed of n transient orbits, and it is labelled by a set of n indices 
71, . . . , 7 n . It is customary to consider cyclic permutations of these indices to correspond to the 
same diffractive orbit. The contribution of a diffractive orbit to the density of states in (|T9|) can 
be interpreted in the following way. In the GTD approximation the diffraction is treated as a 
local process that occurs at the source of diffraction. The diffraction coefficient 2?(0 , 0& >7i ) 
contains the amplitude and phase for the scattering from the incoming direction of the transient 
orbits 7j into the outgoing direction of the transient orbit 7^+1. In some sense, the diffractive 
trajectories can be considered as generalised periodic orbits which are closed in momentum 
space by the scattering at the singularity of the potential. 

Some examples for the diffraction coefficients are 



Diffraction on a corner with angle 6 in a two-dimensional billiard system (ITO, IT? 
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where N = 9/ir, and (f> a and </>(, are the angles of the outgoing and incoming trajectory, 
measured with respect to one side of the corner such that (fi a , 4>b £ [0, 0]. 



Diffraction on a flux line in a two-dimensional system [I8|: 
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where a is the flux parameter and <p a and <pb are the directions of the outgoing and 
incoming trajectory. 



Diffraction on a delta-like singularity of the potential in two dimensions flpj, EDI : 
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where k = \/2mE/% and a is a parameter characterising the strength of the potential, 
and 7 is Euler's constant. 



Diffraction on a delta-like singularity of the potential in three dimensions [IT5, 2U|: 
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In the first two examples the diffraction coefficient is energy independent but it depends 
on the incoming and outgoing directions of the trajectories. There are directions in which 
the diffraction coefficient diverges. There the GTD approximation breaks down and has to be 
replaced by a uniform approximation |2l], |l8fl . In the last two cases, the diffraction coefficient 
has no angular dependence, but depends on energy. It corresponds to pure s-wave scattering, 
and the GTD approximation is valid for all angular directions. 

In the following we consider the influence of diffractive orbits on semiclassical approxima- 
tions for the spectral form factor which is defined by 
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where (■ ■ -)e denotes an average over an energy interval that is small in comparison with E but 
contains many energy levels. If the oscillatory part of the density of states is semiclassically 
approximated by classical trajectories in the form 
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then the spectral form factor is expressed by a double sum over trajectories. We consider here 
only the diagonal approximation to this double sum that describes the form factor for small 
values of r. In leading semiclassical order 
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where T = 2irhd(E)T and g 1 is the number of terms in the sum for which the actions S y (E) 
are identical. A 7 is the average over all g 1 amplitudes of the orbits which have the same action 
as the orbit 7 in case that these amplitudes are different. 

Let us first consider the case of single- diffractive orbits n = 1. Then 
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and one obtains for the contribution of diffractive orbits with one scattering event to the 
diagonal approximation of the form factor 
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The degeneracy of the actions is in this case g = 2/(3 (except for an exponentially small fraction 
of the orbits for large T) where (3 is an integer denoting the symmetry class of the system. (3 = 1 
for systems in which the only symmetry is an anti-unitary symmetry, and (3 = 2 for systems 
without any symmetry. Applying the sum rule flT5| ) and the leading asymptotic approximation 
for the mean density of states d(E) ~ Yi(E) / (2itK)f one obtains 
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t) from diffractive orbits with fixed initial and final direc- 



for the partial contribution to K d 
tions. 

For the examples given above this result has the following implications. For the corner 
diffraction and diffraction on a flux line the prefactor of r 2 is independent of h. It implies 
that diffraction has an influence on statistical properties of energy levels in the semiclassical 
limit h — *■ 0. Although the contributions of diffractive orbits to the density of states are by 
an order \fh smaller than those of periodic orbits, they still give a finite contribution to the 
diagonal approximation of the spectral form factor as h — > 0. The complete contribution to 
K d that is obtained by integrating over all angular directions <p a and <fi b requires the use 
of uniform approximations, since the GTD approximation yields a divergent result. However, 
since all contributions to K d are positive there is a non- vanishing contribution of order r 2 to 
the diagonal approximation of the form factor in the semiclassical limit. For the examples of 
diffraction on delta-singularities the diffraction is isotropic and the integration over the angular 
directions can be performed as is done below. Here the prefactor of r 2 is energy dependent, 
and there is a significant difference between the two- and three-dimensional result. In two 
dimensions the prefactor vanishes in the semiclassical limit, although very slowly like (log ft) -2 , 
whereas in three dimensions the prefactor approaches a constant. 

The contributions of diffractive orbits with n scattering events to K<i(t) are obtained in 
a similar way. They are composed of n transient orbits and it is sufficient to consider only 
cases in which all these n orbits are different, since the relative number of the other cases is 
exponentially suppressed for large times T. The degeneracy of the actions of diffractive orbits 
is then g = (2/(3) n n\. The index 7 in (B^) is now a multiple index 7 = (71, . . . , 7 n ) and 
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is the average over all n\ permutations of the 7$. For simplicity of notation we abbreviate the 
argument of T> by an index 7. Using the relation 
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one can write the contribution of n-fold diffractive orbits to the diagonal approximation of the 
form factor in the form 



Kf\r) 



TH2/P)"n\ r° ^ 



i=\ 



det' 



dq b 
dp a 



m - T 7 



(47rmv a ) 2 (2irHy- 3 
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which can be shown by induction, and the final result is 
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The interpretation is similar as before. The n-fold diffractive orbits contribute to Kd{r) in order 
t u+1 . For corner diffraction and diffraction on a flux line the prefactor of r n+1 is independent 
of h and persists in the semiclassical limit. The remarkable point is that the order in ft of the 
contributions of these orbits to the density of states can be arbitrarily large, their amplitude is 
by an order h n l 2 smaller than those of periodic orbits, but they still give a finite contribution 
to the form factor in the semiclassical limit! For the diffraction on delta- singularities of the po- 
tential, the prefactor of r n+1 again vanishes slowly in the two-dimensional case and approaches 
a constant in three dimensions as h — > 0. 

For the isotropic case where the diffraction coefficient is angular independent, the integration 
over the angular directions can be performed, and the contributions of all diffractive orbits can 
be summed. Adding also the contribution of periodic orbits, the complete expression for the 
diagonal form factor is 
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For an s-wave scatterer in a three-dimensional billiard system the constant approaches the value 
C = 8 in the limit h — > 0. 

The results for the contributions of diffractive orbits to the spectral form factor can also be 
applied to a rectangular billiard with an s-wave scatterer, since the same sum-rule applies in 



this case. The form factor for this system has been studied in |]22 
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4 Conclusions 

The main result of this article is that diffraction can have an influence on spectral statistics 
in the semiclassical limit. The semiclassical treatment of diffraction leads to corrections to the 
diagonal approximation for the spectral form factor K(t) that do not vanish in the semiclassical 
limit h — > 0. This is in contrast, for example, to corrections due to bouncing ball orbits, which 
vanish if the semiclassical limit is performed while keeping the argument of K(t) fixed. 

Although the corrections to the form factor do not agree with those expected from random 
matrix theory, the present results do not show that diffraction leads to a deviation from random 
matrix statistics, since off-diagonal contributions to the form factor have been neglected. They 
show, however, that a semiclassical proof of the random matrix hypothesis cannot be based on 
periodic orbits alone if diffraction occurs. Consider, for example, chaotic billiard systems with 
flux lines. These systems are considered to be standard examples in which spectral correlations 
are expected to follow the statistics of the Gaussian unitary ensemble (GUE). The diagonal 
approximation for the form factor in terms of periodic orbits already yields the correct linear 
GUE form factor up to the Heisenberg time r = 1 |. Moreover, it can be shown that in 



certain ensemble averages the off-diagonal contributions of periodic orbits vanish [23|. The 
present results show, however, that diffraction on the flux line leads to additional semiclassical 
contributions to the diagonal form factor in all higher orders of r. If the spectral statistics 
follow indeed random matrix theory, then these additional contributions have to be cancelled 
by off-diagonal terms involving diffractive orbits. Such a cancellation would require specific 
correlations between classical trajectories that start from and return to one point in coordinate 
space, or between these trajectories and periodic orbits, analogous to the action correlations of 
periodic orbits J7|. 

For another example let us assume that the spectral correlations in a particular three- 
dimensional chaotic system are described correctly by random matrix theory in the semiclassical 
limit. Then it is generally expected that this property is not changed, if an s-wave scatterer is 
added to the system. As is shown in the present article, however, the diagonal approximation 
for K(t) is modified also in this case. This would imply a deviation from random matrix results, 
if this is not corrected by the off-diagonal terms. Nevertheless, it shows that a single s-wave 
scatterer can lead to deviations from random matrix statistics. 

The analysis in this article is based on the geometrical theory of diffraction. For corner 
diffraction or diffraction on a flux line this is not sufficient for calculating the complete diagonal 
approximation of the form factor. It would require the application of uniform approximations. 
Since the semiclassical weight of diffractive orbits is larger in the uniform regime, one can expect 
an even stronger total influence of diffraction on the spectral form factor in these cases. The 
main remaining question is, whether deviations from random matrix statistics can be observed 
in these systems. 
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A Cat maps 

In this appendix we give a simple example for the sum rule for transient orbits. Cat maps are 
linear, area-preserving, hyperbolic maps on the unit 2-torus 

^ n+1 ) = M (f) mod 1 , (39) 

where detM = 1 due to area preservation, | TrM| > 2 due to hyperbolicity, and the elements 
Mij of the matrix M are integers for continuity. 
The sum rule (|i~3"D has the form 

, £ ,ra =1 ' (40) 

where the sum extends over all points (qo,Po) for which q = q a and q n = qv The matrix 
element (M ra ) 12 is the same for all points and can be taken in front of the sum. The number of 
points over which the sum extends is given by the number of solutions of the equation 

q b = (M n ) u q a + (M n ) 12Pa modi. (41) 

Since p a varies in the interval [0, 1) there are exactly |(M n )i2| solutions of this equation. This 
shows that the sum rule ([40| ) is exact for all n > 0. 

Let |TrM| = 2cosh(A) where A > is the Lyapunov exponent of the map, and o = 
sign(TrM). Then 

M* = I h S A f^ enX ° \( s * ~ S A (42) 

~ u 1 s 2 -u 2 s 1 \u2 s 2 )\ a n e~ nX ) \-u 2 Ul J ' l } 

where (ui,u 2 ) and {si,s 2 ) are the components of the unstable and stable eigenvectors of M, 
respectively. One obtains (M n ) 12 = — <r n UiSi(e nX + e~ nX ) / \u\s 2 — u 2 S\) and it follows that the 
number of points contributing to the sum rule is given by 

p(n) ~ - — t^il — - e nX , n -» oo . (43) 

\U\S 2 — U 2 Si\ 



This is of the form (18) and it shows, moreover, that the constant c in fll8|) is not universal 



For cat maps it is determined by the directions of the eigenvectors of M. 
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